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Abstract 

Two-loop contributions to the electromagnetic form factors are calculated in the kinematic 
regime close to the fermion-antifermion threshold. The results are presented in an expansion 
in the velocity /3 of the fermions in the cm. frame up to next-to-next-to leading order in p. 
The existence of a new Coulomb singularity logarithmic in (3, which is closely related to the 
0{a 2 In a) corrections known from positronium decays, is demonstrated. It is shown that due 
to this Coulomb singularity 0(a 2 ) relativistic corrections to the non-relativistic cross section 
of heavy fermion-antifermion pair production in e + e~ annihilation cannot be determined by 
means of conventional multi-loop perturbation theory. 



1 Introduction 



In view of future experiments (NLC, B-factory, r-charm factory) where heavy quark-antiquark pairs 
will be produced in the kinematic region close to the threshold and a large amount of data can be 
expected, it is a very attractive idea that an extraction of the strong coupling a s at a specific scale 
(or equivalently A QCD ) might be possible which is accurate enough to allow for a serious comparison 
to complementary determinations of a s from high energy experiments, where quark masses are much 
smaller than the relevant energy scales. Such an analysis would be an extremely important test of 
QCD. In recent literature two attempts can be found [M ^] where such an analysis has been carried out 
based on present data on properties of bb mesons and on theoretical calculations involving well known 
results in the non-relativistic limit. The results of these analyses are somewhat controversial indicating 
that a better understanding of the structure and size of relativistic corrections to the non-relativistic 
limit and of the interplay of these corrections with non-perturbative effects is mandatory. 

The framework in which relativistic corrections can be determined systematically in a very 
elegant way is non-relativistic quantum chromodynamics (NRQCD) |J which is based on the concept of 
effective field theories. NRQCD consists of a non-relativistic Schrddinger field theory with a Coulomb- 
like QCD potential whereby relativistic effects are incorporated by introduction of higher dimensional 
operators in accordance to the underlying symmetries. In order to render NRQCD equivalent to QCD 
the NRQCD Lagrangian has to be matched to predictions in the framework of conventional multi- 
loop perturbation theory. This procedure leads to, in general, divergent renormalization constants 
multiplying the operators in the NRQCD Lagrangian and is essentially equivalent to a separation of 
short- and long-distance effects. As far as the decay and production properties of a heavy quark- 
antiquark pair involving single photon annihilation in the threshold regime are concerned the relevant 
parts of the NRQCD Lagrangian have only been renormalized at leading and next-to-leading order in 
a s so far 0. 

In this letter we present the two-loop contributions to the electromagnetic vertex describing 
the decay of a virtual photon into two massive fermions in the kinematic regime where the squared 
photon four momentum is close to four times the squared fermion mass. The calculation is performed 
in the framework of QED where only one fermion species with mass M and electronic charge e exists. 
The result is presented up to next-to-next-to-leading order (NNLO) in an expansion in 

/ M 2 

f3 = Jl-A^—, (1) 
y q z + ie 

which is equal to the velocity of the fermions in the cm. frame above threshold]], \fq 1 ~ being the cm. 
energy. We analyse the structure and form of the results and demonstrate the existence of a new 
logarithmic Coulomb singularity occurring at NNLO in the velocity expansion. In particular, we will 
study the impact of this singularity on the massive fermion-antifermion pair production cross section 
slightly above the threshold. In the framework of QCD our two-loop results represent all two-loop 
contributions involving the color factor Cp (from exchange of two virtual gluons) and CfT (from 
the exchange of one gluon with the insertion of the fermion-antifermion vacuum polarization) and, 
therefore, are a gauge invariant subset of all two-loop QCD contributions in the threshold regimq^j. 

1 Thus f3 will be called "velocity" for the rest of this paper. In this paper we use the notion "leading order" (and 
NLO, NNLO, NNNLO) exclusively for the expansion in the velocity. 

2 The two-loop contributions arising from the virtual effects of massless fermions have been calculated in H for all 
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The two-loop contributions calculated in this work represent a first step toward a two-loop 
renormalization of the NRQCD Lagrangian describing single photon annihilation processes involving 
heavy quark-antiquark pairs. In particular, they are a crucial input for the determination of NNLO 
relativistic corrections for the single photon annihilation contributions to decay and production of 
heavy quark-antiquark bound states and for the production of heavy quark-antiquark pairs in e + e~ 
collisions slightly above threshold. In the framework of QED the result is essential for the determina- 
tion of the single photon annihilation contributions to the (D(a 6 ) triplet-singlet hyperfine splitting of 
the positronium ground state. 

The content of this work is organized as follows: in Section |2| we explain the notation and 
introduce the electromagnetic form factors relevant for our calculations and discussions. In Section ^ 
we reanalyse the well-known one-loop contributions to the form factors in the threshold region. We 
discuss the structure and properties of the individual coefficients of the expansion in small j3 and 
derive predictions for the form of the two-loop corrections based on the factorization of long- and 
short-distance contributions. In Section [| the two-loop corrections are explicitly calculated using the 
dispersion integration technique. It is demonstrated that the predictions of Section |3] are realized and 
the logarithmic Coulomb singularity is discussed. Section || contains a summary. 



2 Notation and Definition of the Electromagnetic Form Factors 

It is common to parameterize radiative (multi-loop) corrections to the electromagnetic vertex, describ- 
ing the decay of a photon with virtuality q 2 into a fermion-antifermion pair, in terms of the Dirac (F\) 
and the Pauli (F2) form factors. They are defined through the relation 



u(p') h e J n v(p) = ieu(p') 



J f .F 1 (q 2 ) + —a flu q»F 2 (q 



v(p) , 



(2) 



where 



q = p + p' and 



a 



fJbV 



Expanded in the number of loops, which corresponds to an expansion in powers of the fine structure 
constant a, the form factors F\ and F2 read 



fl(9 2 ) 



i + (f)**V 



+ 



+ 



a\2 



(?) 



a\ z _(2) 

2 



+ 



(3) 



The use of F\ and F2 is particularly convenient for the kinematic point q 2 = because ^(O) = (gf — 
2)/2 is directly related to the gyro-magnetic ratio of the fermion and because i*i(0) = 1 (i.e. F^ n \o) = 
for n = 1,2, ...,00) due to gauge invariance. These properties are useful if dispersion relation 
techniques are used to calculate higher loop contributions because overall UV divergences to F^ 
(n = 1,2, ...,00) can be automatically renormalized by using once-subtracted dispersion relations. 
For F2, on the other hand, no overall UV divergences exist which makes the use of unsubtracted 



ratios M 2 /q 2 above threshold and will not be discussed in this work. 



3 



dispersion relations convenient. Since the determination of our two-loop results relies on the dispersion 
relation technique we will use the form factors F\ and F% for the actual calculations. 

For physical applications in the threshold region, where q 2 ~ 4M 2 , however, the use of the 
combinations 



G m 
G e 



F x +F 2 



F x + 



AM 2 



F, 



(4) 
(5) 



is more appropriate. This can be easily seen by considering the contributions of the form factors F\ 
and F% to the cross section for the production of a fermion-antifermion pair (with fermion mass M) in 
e + e~ annihilation. Taking the colliding electrons and positrons as massless one arrives at the following 
angular distribution for the produced fermion pairs for the cm. energy \fq 2 above threshold 



dcr(e^ 



dO 



ff) _ a^_p 
~ Tq 2 



\G r . 



Jl + cos 2 6») + 



4M 2 



\GA 2 sin 2 1 



(6) 



where 9 is the deflection angle. The corresponding expression for the total cross section reads (a p t 
4ira 2 /3q 2 ) 



R = 



aye' e 



ff) 



P 



\G r , 



P 2 ) \G e 



(7) 



G m and G e are called magnetic and electric form factors, respectively ||. They can be easily identified 
as the total spin projection (relative to the electron direction) ±1 and amplitudes describing the 
produced fermion-antifermion pair in a triplet (J = 1 ) state. Because the fermion-antifermion 
production cross section represents one of the most important applications of the corrections to the 
electromagnetic vertex we will discuss the structure and properties of the corrections by analysing the 
moduli squared of the magnetic and electric form factors above threshold. Their expansion in the 
number of loops (i.e. in powers of the fine structure constant) reads 



1 a, 



\G, 



1 + 



7T 



+ 



+ 



-^ 2 e + 



TV' 

a\2 



(J, 



(2) 



+ 



(8) 



We finally would like to emphasize that throughout this paper the fermions are understood as stable 
particles and that the on-shell renormalization scheme is employed, where a = 1/137 and M is the 
fermion pole mass. 



3 One-Loop Results 

Analytic expressions for the one-loop contributions to the electromagnetic vertex valid for all energies 
are well known since quite a long time @, |j. In this section we reanalyse the one- loop contributions 
in the threshold region in the velocity expansion as a preparation for the examination of the two-loop 
contributions in Section ||. 



4 



Regularizing the soft photon infrared divergences with a fictitious small photon mass A, where 
the hierarchy X/M <C \(3\ <C 1 is understood, the one-loop contributions to the electromagnetic form 
factors F\ and F2 assume the form 



IT 

4 
3 



ln( 



2ij3M 
A 

M\ 5 
Vt) + 24 



3 7T/3 

- + I — 

2 2 



ln( 



2i(3M 
A 



In (— ) 



■ 7T 1 . vr /? 1 2 
2 47" 2 _, T- + 3* 



/3 2 + 0(/? 3 ) , 
0((3 3 ) 



(9) 
(10) 



in the velocity expansion up to NNNLO. Expressions (||) and (|T^) are valid above as well as below the 
threshold point, q 2 = AM 2 , and lead to the following one-loop contributions to the moduli squared of 
the magnetic and electric form factors above the threshold 



/3^o air 
V3 



a 



TT 



air (3 
2 

an (3 



{-)9 e {q ) = - 4 - + — — 



a 


8 In 


3n 




8a 


ln( 




3tt 







1 " 




~ 3. 



+ 



(3 2 + 0(/3 3 ) 
p 12 + 0((3 3 ) . 



(11) 
(12) 



For the rest of this section we will discuss the individual terms in the velocity expansion displayed in 
eqs. (|Tl| ) and (|l^). We would like to emphasize that most of the issues which are mentioned are well 
known and have been noted before at various places throughout the literature. However, we think 
that a review of these topics is necessary for a better understanding of the structure of the two-loop 
results presented in Section || and the new information contained in them. 

Expressions ( |Tl"l ) and ( |i"2| ) exhibit the well known soft photon divergence oc ln(M/A) which 
arises from the masslessness of the photon. This divergence occurs at order 1 and would cancel with 
the corresponding soft photon divergence coming from the process of real radiation of one photon 
off one of the fermions according to the Kinoshita-Lee-Nauenberg theorem 0, |l(]] . The fact that the 
divergent term ln(M/A) is suppressed by /3 3 relative to the leading contribution in the expansion in (3 
is expected at any loop level because close to threshold the real radiation of one photon results in an 
additional factor j3 from the phase space needed for the photon and a factor (3 2 from the square of the 
dipole matrix element^. Because the soft photon ln(M/ A) divergence indicates the inadequacy of a pure 
fermion-antifermion final state and the need for the introduction of a higher fock fermion-antifermion- 
photon state, the /? 3 suppression allows us to conclude that the notion of a pure fermion-antifermion 
state is consistent if we are only interested in NNLO accuracy in the expansion in (3. 

The leading term in the velocity expansion in eqs. (11) and (|l^) is the well known Coulomb sin- 
gularity which diverges for /3 — ► 0. Similar to the soft photon divergence discussed above the Coulomb 
singularity arises from the fact that the photon is massless and represents a long-distance effect. The 
Coulomb singularity, however, is of completely different nature. Whereas the soft photon singularity 
indicates the inadequacy of a pure fermion-antifermion state beyond NNLO in the velocity expansion 

3 It should be noted that this statement is equivalent to the fact that contributions from the non-instantaneous (i.e. 
transverse) exchange of photons among the fermion-antifermion pair are suppressed by f3 3 with respect to the leading 
contributions in the velocity expansion. As an example, this feature is apparent in a 3 Si, J PC = 1 fermion-antifermion 
bound state, where the velocity /3 of the fermions is of order a. There, the exchange of non-instantaneous photons leads 
to the Lamb shift which represents a 0(a 3 ) correction relative to the Coulomb energy levels. (See also fl~L|.) 
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the Coulomb singularity reveals that in the non-relativistic limit (corresponding to the leading order in 
the velocity expansion) the photon-mediated interaction between the fermion-antifermion pair cannot 
be described in an expansion in Feynman diagrams, where a diagram with a larger number of loops 
(corresponding to a larger number of exchanged photons) would represent a higher order correction. 
Rather, a resummation of diagrams with any number of exchanged photons is needed to arrive at a 
sensible description of the interaction between the fermion-antifermion pair. The leading contribu- 
tion in the velocity expansion is obtained by resummation of diagrams with instantaneous Coulomb 
exchanges of longitudinal photons (in the Coulomb gauge). This procedure can be explicitly carried 
out by calculating the normalized wave function at the origin, ^^O), to the Schrodinger equation 
describing a non-relativistic fermion-antifermion pair with a Coulomb interaction potential for positive 



energies E = M0 2 . The result of this calculation reads (see e.g. f§ 0, Q) 

l G m|LO = I^eIlo = l*M/3 2 (°)| 2 = i _ OVT ^ r~7i ' ( 13 ) 

where 

z 55 T' (14) 

and is often called "Sommerfeld factor" in the literature. The 1/(3 Coulomb singularity in eqs. (11) 



and (12) can be recovered as the 0(a) contribution in the expansion of the Sommerfeld factor for 



a <C (3, 



z 2 



l-exp(-z) 1 + 2 + T2 + °^- ™ 

This, on the other hand, also shows that the velocity expansion of the perturbative (in the number of 
loops) series can only be applied in the limit a</}<l, where an expansion in the number of loops 
(i.e. in a) is justified^. It is worth to study the effect of this resummation: inserting the Sommerfeld 
factor into the formula for the cross section, eq. ([?]), we get at threshold 

d 3 a z 3 na\ 

K ~ -p - — > -air, (lb) 

2 H 1 - exp (-z) 2 v ; 

which is the correct result according to non-relativistic quantum mechanics. On the other hand, if we 
naively use the one-loop result (i.e. expansion in small a), we obtain 

!/?(! + £) (17) 

Clearly, the perturbative calculation in the number of loops, which is based on the assumption that a 
is a valid expansion parameter close to threshold, gives a prediction for R at threshold which deviates 
from the correct one by a factor of one half. 

The next-to leading contribution in the velocity expansion in eqs. (|IT| ) and (|l^), — 4q/7t, 
represents a short-distance correction and can be understood as a finite 0(a) renormalization of 



4 It should be noted that the region of convergence of the Taylor expansion 

oo 

>+?+»-> 



l-exp(-z) * ' 2 L > (2n)l 



where B n are the Bernoulli numbers (Si = 1/6, B 2 = 1/30, B 3 = 1/42,...), is \z\ < 2ty <^> |/3| > a/2. This shows 
that for phenomenological applications a resummation of the leading order contributions in the velocity expansion to 
any number of loops is mandatory in the kinematic regime |/3| < a. 
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the electromagnetic current which produces the fermion-antifermion pair in the threshold region. The 
short-distance character of this 0(a) correction has been demonstrated explicitly by the calculation of 
the BLM (Brodsky-Lepage-Mackenzie [14|]) scale in the coupling governing the — Aa/n contribution 
|l~5| , [lj. This BLM scale is of order the fermion mass M and indicates that the — A a /it contribution 
represents a correction to the fermion-antifermion production process which occurs at short distances 
of order 1/M. In contrast, the BLM scale of the coupling in the leading term in the velocity expansion, 
air /2f3, is of order of the relative momentum of the fermion-antifermion pair, Mj3 p, |15| , indicating 
that the latter contribution belongs to the fermion-antifermion wave function. As a consequence the 
leading order (long-distance) contributions contained in the Sommerfeld factor and the short-distance 
corrections are expected to factorize which leads to 



O' 

1 — exp (—z) V vr 



|G>)iInlo — l^elNLo — i 7 IX ( 1 — 4— ) (18) 



for the NLO expressions in the velocity expansion of the moduli squared of the magnetic and electric 
form factors in the threshold region. It should be noted that the factorized result (|l8|) resums all con- 
tributions (a/(3) n x [1, a], n = 0, 1, 2, . . . , oo. Because no {a/ (5) n f5 contributions existQ, expression (18) 



unambiguously predicts the leading and next-to-leading order contributions in the velocity expansion 
for all g£J e , n = 2, 3, . . . , oo. 

The NNLO term in the velocity expansion in eqs. ( |TT1) and (|i~2|), airj3/2, has not received much 
attention in the literature so far. Its structure, which involved the same power of tt and the same 
coefficient 1/2 as the LO term in the velocity expansion, strongly implies that it is of long-distance 
origin and therefore belongs to the Sommerfeld factor. This is in accordance to the observation that 
the BLM scale in the coupling of the term an (3/ 2 is of order M/3 rather than M ||. The relativistic 
extension of the Sommerfeld factor (including 0([3 2 ) corrections) should then read 

^(l + /3 2 )- (19) 



1 — exp(— z) [5 

Although the arguments given above in favor of expression (|l9| ) are far from being a strict proof 
the form of z is very convincing because it indicates that the relativistic relative velocity v re i of the 
fermion-antifermion pair in the cm. frame is involved in the argument of the Sommerfeld factor if 
0{(5 2 ) relativistic corrections are taken into account, 

2avr 2(5 , . 

Combining expression ( |l9| ) with the short-distance factor (1 — 4a /it) and taking into account that no 
soft photon divergence oc ln(M/A) arises up to NNLO in the velocity expansion we can now predict 
that the two- loop contributions to \G m / e \ 2 must have the form 



4 2 4 

(2) , 2\ /3^0 IT TT TT 



finite terms without 7r 4 



+ O(P) . (21) 



We want to emphasize the the 0(1//? ), C(l//3) and 0(f3 it ) contributions on the r.h.s. of eq. (21) 



are an unambiguous prediction and have to be recovered in the explicit two-loop result if the concept 

5 Pure /3-dependent corrections to the Sommerfeld factor are of kinematic origin and therefore expected to be of 
NNLO in the velocity expansion, i.e. oc (a//3) n /3 2 , n — 0,1,2, ... , oo. 
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of factorization in the threshold regime is valid. It should be noted that up to NNLO in the velocity 
expansion only the 0{fl ) contributions symbolized by [finite terms without 7r 4 ] contain new two-loop 
information. 



4 Two-Loop Results 



To determine the two-loop contributions to the electromagnetic form factors F\ and F2 in the ve- 
locity expansion we use the dispersion integration technique. For that we have to integrate over 

(2) 

the absorptive parts ImFy^ which have been determined a long time ago by Barbieri, Mignaco and 
Remiddi [16], 



ifV) 



2 „2 



AM z q 



p 

AM 2 1 



'(2) 



(0) 



+ ■ 



7T q z 



AM 2 



4M 2 



dq 



12 



,'2 



(q>< 



ie 



-/2 



AM 1 



V2 



lmFl'(q 



(22) 



ifV) 



AM A 



AM 2 F2 



(2) 



(0) 



1 

+ - 



7r q* 



AM 2 



4M 2 



■J 2 



J 2 



AM 2 



-.12 



IE 



,12 



lmF^(q' 2 ) 



(23) 



We would like to mention that relations (|2^) and (^) are equivalent to the common once-subtracted 
and unsubtracted dispersion relations. We use (|22|) and (p3|) because they do not run into non- 
analyticity problems in the integration region where q' 2 — AM 2 is of order A 2 if the limit A — > 
is already taken before the integration. Since the absorptive parts in [|i]] are given in exactly this 



limit (|2^ ) and (23) are more convenient because in them the integration regime q' 2 — AM 2 ~ A 2 is 



strongly suppressed. The (low) price one has to pay is that the 0(a 2 ) fermion charge radius [O, Hi 



F{ (2) (0) 



1 

M 2 



6 



3 In 2 



49 
72 



4819 
5184 



(24) 



and the 0(a 2 ) anomalous magnetic moment [19. 

■2 



if } (0) 



7T" 

12 



3 197 
fi,ll2 - : - 1 l + 4 C3+ 144 



(25) 



have to be taken as an input^j. Details for the quite lengthy but straightforward calculation of the 
integrals ( [22] ) and (p3|), which requires strong support of algebraic manipulation programs, shall be 
presented elsewhere. 

6 This fact has already been pointed out in [[L6|. We also refer the reader to this reference for a more thorough 
discussion of the problems which occur in the integration region q — AM 2 ~ A 2 . 
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The final results for the two-loop contributions to F\ and Fi up to NNLO in the velocity 
expansion read 



F 



(2) ^0 



7T _ 



1,27 



i/3 2 



7T 



+ 



+ 



+ i 



4/3 



3£ + l 



7T 7T 



24 4 1 



23 , , „. 7 , 73 \ 9 / 421 \ 

+ — ln(-i/3) + — ln2 + — + — [9C 3 

15 v HJ 10 50 J 80 V 27 J 



,(2) /?->0 13 7T 2 37 



1,/ 



7 (2) 
2,2 7 - 



+ — + 0{[5 2 ) 
45 12 v ; 



(26) 
+ 0(/3) , 

(27) 



7T 



2 r 



i/3 2 



7T 



4/3 L 



+ 1 



+ 



ttV, . .„ 101 . n 559 \ 1 / , 269 

— ln(-i /3) H In 2 M 41 (3 H 

20 V V H> 6 45 J 80 V S 3 



+ O03) 



,(2) ^0 vr 2 23 



2,/ 



where 



In 



V A ) 



(28) 
(29) 

(30) 



and, as in the one-loop case, the hierarchy X/M <C \/3\ <C 1 is understood. In eqs. fl2q)-(|29D the 
contributions from diagrams with two photons (subscript 27) and from the diagrams with one photon 
and the insertion of the fermion-antifermion vacuum polarizationT] (subscript /) are displayed sepa- 
rately. This will facilitate the application in the framework of QCD where both types of contributions 
are multiplied by the different color factors Cp and CfT, respectively, and represent gauge invariant 
subsets of the full QCD two-loop contributions. 



The results fl2q)-(29) lead to the following two-loop contributions to the moduli squared of the 
magnetic and electric form factors above threshold up to NNLO in the velocity expansion 



71 



7T 2 7T 4 



12/32 2 p + 6 + ^ 



2 „ 4, 29 

-ln/3 + -ln2 

3^3 12 



527 

C3 + ^ + O(0) , (3l) 



7T 



7T 2 7T 4 



12 /3 2 



2 j + y + * 



?ln/3+^ln2-I)-Cs + ^ + 0(«. (32) 



It is evident that the prediction made in the previous section based on the one-loop corrections and 
on the factorization of long- and short-distance contributions (see eq. (^) ) are indeed realized by our 
explicit two-loop result confirming the statements given in Section]^. As a consequence only the 0(f3°) 
terms in eqs. ( |3l"D and (32) essentially contain new information. 

The most conspicuous feature of the O(/3 ) contributions in eqs. ([H]) and ( |32[ ) is the term ln(/3). 
Similar to the 1//3 2 Coulomb singularity exhibited in the leading term in the velocity expansion, it 
indicates the breakdown of the conventional perturbation series in the number of loops in the limit 
(3 — > 0. The existence of this logarithm can be understood from the fact that two scales are involved 
in the kinematic regime near threshold, the fermion mass M and the three momentum of the fermion 
and antifermion in the cm. frame p = M(3. The logarithm of the velocity (3 is therefore actually 
the logarithm of the ratio of these two scales, m(p/M). Because the soft scale p is characteristic 

7 The two-loop corrections F^ 2 J and F^ 2 j have already been calculated before in ffl for all energies above threshold. 
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for the fermion-antifermion wave function and not relevant for the production mechanism of the 
fermion-antifermion pair (which involves only the hard scale M), the a 2 \n(p/M) term in eqs. (31) 
and (^) should occur with the same coefficient in the 0(a 2 ) corrections to the positronium decay 
rates. For a viable comparison, however, we also have to include the fermion-antifermion vacuum 
polarization effects coming from the fact that the fermion-antifermion pair, which is in a J PC = 1-- 
state, can virtually annihilate into one photon. This can be easily achieved by multiplying \G m / e \ 2 
by the factor |1 + XT | 2 , where II is the one-particle-irreducible vacuum polarization function. The 
two- loop contribution to II also contains a logarithm of in the velocity expansion This leads 
to the additional contribution a 2 ln(/3) which has to added to — 2 a 2 ln(/3)/3 from \G m / e \ 2 . (Actually 
the spin average of the logarithmic terms in expressions ( |3~I| ) and (|32| ) has to be taken. This trivially 
results in — 2 a 2 ln(/3)/3 because the logarithmic term is universal in both spin amplitudes.) Because 
the relative momentum of the electron-positron pair in the positronium is of order Ma we can expect 
that the 0(a 2 ) corrections to the ( 3 Si, J PC = 1 ) orthopositronium decay rate should contain 
the contribution a 2 ln(a)/3. This logarithmic 0(a 2 ) correction has indeed been found by explicit 
calculations of higher order correction to the orthopositronium decay rate |22] . We therefore have to 
conclude that the ln(/3) term in eqs. (|3l| ) and (|32|) represents a new type of Coulomb singularity which, 
similar to the power-like 1//3™ singularities, requires a resummation of contributions to all orders in 
the number of loops^. How such a resummation has to be carried out for the ln(/3) term in the vacuum 
polarization has been demonstrated in [21]. 

Finally, we want to discuss the impact of the m(/3) singularity on the cross section of fermion- 
antifermion production very close to threshold, see eqs. @ and (|7]). Because the moduli squared of 
the magnetic and electric form factors are multiplied by the phase space factor f3 one might naively 
think that the ln(/3) singularity is suppressed by (3 and does not affect the cross section for j3 — > 0. At 
this point we have to emphasize that the same would then be true for the short-distance correction, 
— 4a/7r, in the one-loop contribution to |G m / e | 2 because the latter also represents a O([3 ) term in the 
velocity expansion (see eqs. ( |il~l) and (jl2|)). However, the one-loop short-distance correction survives 
for (3 — > 0, see eqs. (|i~6|) and (|l8|) . The resolution of this apparent contradiction comes from the fact 
that due to factorization (see eq. (|l8|)) the one-loop short-distance correction is also contained in the 
0(1/ P) term of the two-loop contribution to \G m / e \ 2 where it multiplies the 0(ot) contribution of 
the expansion of the Sommerfeld factor for small a. This contribution does not vanish in the cross 
section for (3 — > and illustrates the mechanism why the one-loop short-distance correction survives 
in this limit. In order to see that something similar happens to the ln(/3) singularity in the two-loop 
results (|3l]) and (|32|) let us have a closer look on the structure of the one- and two-loop contributions to 
the form factors F\ and i*2- It has been shown by Yennie, Frautschi and Suura [24] that the infrared 
soft photon divergences exponentiate completely. Because real soft photon divergences in |G m / e | 2 
occur only beyond NNLO in the velocity expansion (see Section |||) all soft photon divergences which 
arise up to NNLO in the velocity expansion in eqs. (|2^)-(p9|) can be factorized into a divergent phase 
factor which is known as the Coulomb phase. In the moduli squared of the form factors this phase 
drops out. Since the Coulomb phase has to be considered as an intrinsic property of the fermion- 

8 At this point we would like to mention that the logarithmic Coulomb singularity has also been disussed in [^3| in the 
framework of quarkonia decays. However, it is arg ued in ||] (and also in @) that this singularity (called "logarithmic 
infrared divergence" in Q) would indicate that perturbative QCD could not be applied in the kinematic regime close 
to the threshold. We disagree with this conclusion, because we think that this singularity can be treated by a proper 
resummation of contributions to all orders in the number of loops. 
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antifermion wave function, where the relative momentum 2M/3 is a relevant scale, we can expect that 
the divergent phase factor should involve the logarithm of the ratio 2M/3/X. This feature is indeed 
realized because the sum of Born, one-loop and two-loop contributions to F\ and F2 above threshold 
can be rewritten as 



1 + 



a\2 



F, 



(2) 



,,(2) 



cxp 



a 



+ 1 
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2 \P 

TT 



1-2- T*1J 
1 - 



+ (3 ) e j 
4 ^2 



4/3 



TT 

24 



TT 

20 



23 



ln( 



..^) + Z ln 2 + iiZZ 

M ' 2 90 



+ 



a\2 



TT 

9 

80 



TT 



24 /3 2 
9C 3 -43 



7T _ 



+ 1 
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7T 



7T 



7T 7T 

+ 20 



ln(-*/J) + — In 2 



499 
45" 



+ 



1 

80 



41Cs + 



24 (3 2 
347 



(34) 



The factorized expressions ( |33[ ) and (|34|) predict that at the three-loop level the real parts of the form 
factors Fi and F2 contain the logarithmic and A-independent 0(1/ [3) contributions —23 u 3 tt ln((3) /240/3 
and a 3 7r ln(/3)/80/3, respectively, in the velocity expansion above the threshold. As a consequence, 
|G m | 2 and |G e | 2 both contain the three-loop term — a 3 7rm(/3)/3/3 in the velocity expansion. We 
would like to emphasize that the argument just given cannot be used to determine all three-loop 
contributions, but it clearly shows that a logarithmic Coulomb singularity also exists at order a 3 //3 
which does not vanish in the limit (3 — > in the cross section. The coefficient of this singularity further 
strongly implies that the ln(/3) contributions in |G m | 2 and \G e \ 2 to any number of loops and at NNLO 
in the velocity expansion above threshold can be cast into the factorized form 



\r i 2 



NNLO In /3 — contributions 



1 — exp (—z) 



2 

-a 2 ln/3 

3 H 



(35) 



It is clear from expression (35) and the arguments given above that the logarithmic Coulomb singularity 
does indeed affect the prediction for the cross section for (3 — > 0. In particular, we conclude that 
a conventional fixed order multi-loop calculation is not capable to determine NNLO (i.e. 0(a 2 )) 
relativistic corrections to the non-relativistic cross section^. In order to determine the correct form 
of the NNLO relativistic corrections to the non-relativistic cross section (or the Sommerfeld factor) 
resummations of the type mentioned before have to be performed. Such a program is beyond the 
scope of this work and will be carried out elsewhere. 



9 In a recent publication where large-n QCD sum rules were applied to the bb system 0] it was claimed that 0(a 2 a ) 
accuracy was achieved in the determination of the strong coupling a s and the bottom mass because two-loop corrections 
to the cross section were taken into account. Because the large-n limit peels out the threshold behavior of the bb 
production cross section, the results presented in 0] do not include NNLO relativistic effects properly and, therefore, are 
not at the 0(o? a ) accuracy level. (See also fell.) 
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5 Summary 



In this work we have determined the two-loop contributions to the electromagnetic form factors in 
the kinematic regime close to the fermion-antifermion threshold up to NNLO in an expansion in the 
velocity of the fermions in the cm. frame. In the framework of NRQCD and NRQED the results 
are an important input for the two-loop renormalization of the effective Lagrangian. As the main 
outcome of this work we have demonstrated the existence of a new logarithmic (in the velocity) 
Coulomb singularity at NNLO in the velocity expansion. This logarithmic contribution belongs to the 
fermion-antifermion wave function and exists for the production of free fermion-antifermion pairs above 
threshold as well as for fermion-antifermion pairs in a bound state. For the case of fermion-antifermion 
pair production in e + e~ annihilation the logarithm indicates that a resummation of contributions to 
any number of loops is mandatory in order to arrive at a viable (i.e. finite) prediction for the cross 
section with NNLO accuracy very close to the threshold point. 
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